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equation of state with 3/2 < F < 2 exhibit a curvature singularity in Palatini 
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f{R) gravity as a viable alternative to General Relativity. In the present paper we 
further investigate this characteristic of Palatini f{R) gravity in order to clarify 
its physical interpretation and consequences. 
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1. Introduction 

According to recent cosmological observations jlj, the late time evolution of the 
universe seems to be dominated by a cosmological constant or by some unknown 
form of energy [dark energy) which mimics the behaviour of a cosmological constant. 
The many problems connected with the inclusion of such a constant in the standard 
framework of general-relativistic cosmology [2] have led many authors to consider 
possible alternative explanations for the cosmological data. Clearly, since gravity is 
by far the most important interaction governing the dynamics of the universe on large 
scales, one of these alternatives is to modify the theory of gravity itself by changing 
General Relativity (GR) in some way. 

We focus here on one specific generalization of Einstein's theory: Palatini f{R) 
gravity [3] (see [4] for a recent review of other attempts to generalize GR). As can be 
found in many textbooks (see for example [5]), Einstein's theory can be derived from 
the Einstein-Hilbert action not only by means of the standard metric variation, but 
also by taking independent variations with respect to the metric and the connection. 
In this approach, known as the Palatini variational approach, the metric and the 
connection are treated as independent quantities, and one has to vary the action with 
respect to both of them in order to obtain the field equations. The Riemann tensor 
R^fj^crv ^^'^ the Ricci tensor are defined with respect to the now independent 
connection F^^^^ and do not necessarily coincide with the Ricci and Riemann tensors 
of the metric 5^^. Similarly, the Ricci scalar is defined a,s R — g^'^R^iy- If the 
Lagrangian is linear in R (the Einstein-Hilbert action), variation with respect to the 
independent connection forces it to reduce to the Levi-Civita connection of the metric, 
whereas variation with respect to the metric gives the standard Einstein equations. 
Therefore, in the case of the Einstein-Hilbert action the outcome of Palatini variation 
is standard GR. However, clearly Einstein's theory is no longer recovered for a generic 
action 

where f{R) is a function of i?, g is the determinant of the metric g^jy, Sm is the matter 
action and ip collectively denotes the matter fields. (In this equation, as well as in the 
rest of this paper, we are using units in which c ^ G — I). The resulting theory is 
then what is known as f{R) gravity in the Palatini formalism or simply Palatini f{R) 
gravity. It is easy to see that independent variation of the action H]) with respect to 
the metric and the connection gives 

FiR)R^u-lf{R)9^u^&7rT^,, (2) 
V,[V^F(i?)gn =0, (3) 

where F{R) = df/dR, T^^ = —2{—g)~^/'^5SM/5g^'' is the usual stress-energy tensor 
of the matter and is the covariant derivative with respect to the connection F^^^^. 
Note that a crucial assumption has been made in order to derive © and ([3]): the 
matter action has been taken to be independent of the connection F'^'^^, [see ([T|)]. 
This assumption is physically meaningful because it implies that the connection which 
defines parallel transport, and therefore the covariant derivative of matter fields, is the 
Levi-Civita connection of the metric. This demotes the independent connection to the 
role of an auxiliary field [SUZIIH]- Additionally, under this assumption, the Levi-Civita 
connection becomes the one with respect to which the matter stress-energy tensor is 
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conserved [9^ (which imphes, in particular, that test particles follow geodesies of the 
metric gi^u)- In order to restore the geometrical nature of the independent connection, 
one has to allow it to couple to the matter. This leads to metric-afRne f{R) gravity [5], 
which is a different theory with enriched phenomenology [71 [H] . 

Specific choices for the function f{R) in the action ([!]) have been shown to 
lead to models of Palatini f{R) gravity which might be able to address dark-energy 
problems [10^ . There is now an extensive literature on the cosmological aspects of 
such models and discussing their consistency with cosmological [11] and Solar System 
constraints [TH [T3] . In a previous paper [Tl] we focused on the less well-studied issue 
of finding consistent solutions for static spherically-symmetric matter configurations 
when f{R) ^ R. In order to be able to treat the field equations analytically, we 
assumed a polytropic equation of state (EOS) for the matter, i.e. 

P^kpI, r>3/2 (4) 

[p and po are the pressure and the rest-mass density, while k and T are constants). 
This is a very common and useful choice for making simplified calculations both in 
GR and in Newtonian theory jl5j . We found that for a polytropic index in the range 
3/2 < r < 2 there exist no static and spherically-symmetric regular solutions to the 
field equations, because curvature singularities unavoidably arise at the surface. 
There are four points that ought to be stressed about this result: 

(i) It holds also for any EOS which can be approximated, near to the surface, by a 
polytrope with 3/2 < E < 2. 

(ii) It is independent of the functional form of f{R) (with the exception of some very 
specific cases, including standard GR) [14^. As such, it is applicable not only to 
specific models, but it reveals a generic aspect of Palatini f{R) gravity as a class 
of theories. 

(iii) The singularities appearing are true curvature singularities and not coordinate 
singularities, i.e. the curvature invariants of the metric diverge. 

(iv) Apart from the assumptions already listed, concerning symmetries and EOS, no 
other assumption or approximation has been used. As such, the result applies in 
all regimes ranging from Newtonian weak field to strong gravity. 

As noted in [14] , these results cast some serious doubt on the viability of Palatini 
f{R) gravity. In the next section, after briefly reviewing the arguments justifying 
this claim, we further analyze the situation by considering gedanken experiments 
as a powerful tool to investigate the completeness of the theory (section 12. 1|) . In 
section 12. 2[ we calculate the tidal forces exerted due to the presence of the surface 
singularities, and show that the lengthscale on which they arise is much larger than 
the lengthscale on which the fluid approximation breaks down, unless one considers 
very compact configurations and a very special form for the function f{R) in which 
one cancels by hand several terms generically expected to be present in a cosmological 
scenario. (These two hypotheses were assumed in a restricted version of the calculation 
performed in [16j). In section[31 we then discuss the physical and mathematical nature 
of the problem. This analysis reveals that the presence of the singularities is not 
specifically related to the fiuid description of matter, but rather is a feature of the 
differential structure of the equations of the theory and would, in general, become 
even more acute if the fluid approximation were to be abandoned. In the same section 
we also propose ways to generalize the theory in order to avoid these problems. In 
section |4] we present our conclusions. 
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2. Surface singularities and viability 

2.1. Gedanken experiments and incompleteness 

Clearly, a polytropic EOS is too idealized to give a detailed description for a matter 
configuration resembling an astrophysical star. However, this does not at all make 
polytropes physically irrelevant. On the contrary, as well as being widely used in GR 
and in Newtonian theory [15] . there are at least two physical matter configurations 
which are exactly described by a F — 5/3 polytrope: a monatomic isentropic gas 
and a degenerate non-relativistic electron gas. Note that this value of the polytropic 
index lies well within the range 3/2 < F < 2, for which surface singularities have been 
shown to appear, and so Palatini f{R) gravity does not allow a physical solution for 
these configurations (a solution which is singular at the surface should be discarded as 
unphysical). One might, therefore, discard Palatini f{R) gravity as a viable alternative 
to GR already on the basis of such gedanken experiments. Alternatively, one must at 
least accept that the theory, as it stands, is incomplete, being incapable of describing 
configurations, such as a cloud of monatomic gas, which are well-described even by 
Newtonian gravity. Note that this means in particular that Palatini f{R) gravity does 
not even reproduce the Newtonian limit in these cases! 

It should be stressed that although the fluid description of matter does indeed 
conceal information about the microphysics of the system, this is by no means the 
cause of the problem discussed here, nor will abandoning the fluid approximation 
solve the problem, as we will show in section [3l On the contrary, one naturally 
expects that systems such as a monatomic isentropic gas or a degenerate electron 
gas should be describable by a theory of gravity without resorting to a statistical 
description. In our opinion, the inability of a theory to provide a classical treatment 
of macroscopic systems without a precise microphysical description is already a very 
serious shortcoming. This problem does not arise in standard GR. 

2.2. Stars and tidal forces 

In this section we calculate the tidal forces arising due to the presence of the surface 
singularities which we discovered in [14j . A version of this calculation for a particular 
restricted form of f{R) was performed by Kainulainen et al. [16j . who found that 
the lengthscale on which the tidal forces diverge due to the curvature singularity was 
shorter than the mean free path (MFP) in that case, and concluded that the system 
was not then well-described using the fluid approximation. We will now show that 
while this is correct in the particular case which they considered, that is a very special 
one and is not representative of the general situation. Reference [TB] considered in 
fact the case of a neutron star with f{R) = R — n'^/R (where ~ A, A being 
the value of the cosmological constant as inferred from cosmological observations). 
Although f{R) = R — fi^/R can be used to obtain the accelerated expansion of the 
Universe without resorting to Dark Energy or a cosmological constant, there is no basic 
principle from which to derive this functional form, and in order to justify it one has 
to invoke phenomenological arguments based on a series expansion of the unknown 
f{R) coming from a consistent high energy theory. As such, there is no reason to 
exclude the presence of quadratic or cubic terms, and indeed the observational limits 
on these terms coming from solar system tests are very loose [T7]- We will show 
that if one takes f{R) = R — H^/R + sR^ even with e being orders of magnitude 
smaller than the maximum allowed by the solar system constraints, the lengthscale 
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on which the tidal forces diverge is much larger than the MFP, even in the case of 
neutron stars. Incidentally, this was expected because we have already shown in [T3] 
how important the effect of such a tiny e can be in neutron star interiors. However, 
even if one cancels by hand all of the quadratic and cubic terms from the function 
f{R), thus giving precisely f{R) — R — fi^/R, the result of [16] still does not apply for 
sufRciently diffuse systems, where the lengthscale on which the tidal forces diverge is 
anyway much larger than the MFP. 

Let us first recall the notation and the main results of [14 . In particular, we write 
the static spherically symmetric metric as 

ds^ = -e^(-)dt^ + e^^'-^dr^ + r^dQ^, (5) 

and denote the pressure, energy density and stress energy-tensor of the fluid by p, p 
and T'"'. Also, we define F{R) = df/dR and C = dF/dp{p + p), and use a "prime" 
to denote derivatives with respect to the radial coordinate r. We recall that taking 
the trace of one gets 

F{R)R~2f{R) ^SnT, (6) 

which is an algebraic equation relating R and T for a given /(i?) ^18J . In particular, 
equation ^ implies that in the exterior R is constant [i.e., R = Rq with F{Ro)Rq — 
2/(i?o) = 0]. Similarly, we denote the exterior values of / and F with a "zero": 
/o = /(i?o) and = FiRo). 

Solving ([3]) for the connection and inserting the resulting expression into 12]), it 
is possible to rewrite the field equations in a more familiar form: 

- 1^2 ((V,F)(V.F) - lg,A^Ff^ , (7) 

where and G^^i/ are the covariant derivative and Einstein tensor built with the 
Levi-Civita connection of 5^,^, □ = ^^''V^Viy, and F, f and R are expressed as 
functions of T using ([6]). Because ([7]) reduces in vacuum basically to GR with a 
cosmological constant, it is not surprising that inserting the ansatz ([5]) into it gives 
that the exterior solution must be the well-known Schwarzschild-DeSitter metric {i.e., 
Birkhoff's theorem holds also in Palatini f{R) gravity). This solution is characterized 
by two parameters: the total mass mtot, which is fixed by matching with the interior 
solution imposing continuity at the surface, and the cosmological constant, which 
turns out to be A = Rq/4:. Solving ([7]) in the interior, it can be shown [14] that the 
radial derivative F' is zero at the surface (r = rout), while 

r [rout) - -z 7 7i T ■ i^j 

8rout (rout - 2TOtot ) 

and 

, , , _ 2foflorg,t + {rl^.Ro ~ 8mtot) C 

'^totl^'outj — I6F0 ' 

For 3/2 < r < 2, it can be checked that C ^ 00 as the surface is approached, thus 
driving to infinity F" , m(ot and, more importantly, the Riemann tensor of the metric, 
RfivaX, and curvature invariants, such as R or Ri^""^^ R^^^x. 

We will now proceed to calculate in detail the tidal force experienced, because 
of this curvature singularity, by a body falling radially into a polytropic sphere 
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with r — 5/3, as soon as it crosses the surface. Our conchisions apply unchanged 
also to bodies moving on different orbits, e.g. circular ones just below the surface. 
Throughout the calculation, we use units in which Mq = 1, as well as G = c = 1. 
If we consider the separation vector rj = rfd/dr, the tidal acceleration in the radial 
direction is given by the geodesic deviation equation: 



e^-^(A'^ - A'B' + 2A")(uyr] 



(10) 



where r is the proper time and D / Dt is the total covariant derivative with respect to 
it. Using (8)-(ll) of [T3| and Mathematica [Hj, it is easy to show that the combination 
A'^ — A'B' + 2 A" appearing in this equation depends linearly on F": 

A'^ - A'B' + 2A" = Co + ciF" 



(11) 



where 



Co = {l6F 



+ AOrF'F^ 



[(/ + 16np - FRy + 2F] [p + p)F^ 
2Ar^F'^F ~ 2e^ {2F + rF') [3(/ + 127rp - ATTp)F'r^ 



F (r (r (/' + 2RF') + Stt (p + 
F^ (i?V3 + 4)]F + 3r*F'4|/[: 
4 



p + 2rp')) 
r'^F^ (2F4 



- AF') r 
rF'f] 



(12) 

= -2FTVF ■ 

Note that both co and ci are finite at the surface, whereas F" diverges, as already 
mentioned. Keeping therefore only the divergent term ciF" , the ratio between the 
accelerations in the Palatini ( "singular" ) and GR cases is 



\ciF"\rl^,{r,ut~2mtot) 



(14) 



acR 8mtot 
Using now the fact that Rq must be <C 1 in our units in order to match the cosmological 
accelerated expansion (one needs to have Rq — 4:A — 12QaHq ~ 10"^^), ([8]) gives 
F" w -nitotC /[rout{rout - 2mtot)] and therefore 

|ciC'|rout 



(15) 



The derivative of C with respect to r can easily be calculated from the definition 
C = dF /dp (p + p): using the chain rule, the Euler equation p' = —A'{p + p)/2 and 
the fact that the trace of the stress energy tensor for a perfect fluid is T = 3p — p, one 
has 



C 



dC 
dp 



P 




dF dR 
dRdf 
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Remembering now that C = 0, dp/dp{p + p) ^ and dp/ dp = at the surface for 
r < 2 |14j . one can easily rewrite the above equation keeping only the divergent terms: 

2 



_ dF dR A' 
^ dRdrY 



terms going to zero at the surface 



Taking F = 5/3 and using A' 
equation becomes 



2mtot/[''out(rout - 2mtot)] (see (17) of [14^) 



(17) 
this 



C 



dF dR A' 
dRdrY 

dF 
-Sn- 



ip +P) 



dp 
dp 

3rntot 



dp'^ 



P? 



)/t2 (k) 



-1/5 



di? 25rout(r-out - 2mtot)'"'^ v .. / ' ( ^) 

where, in order to pass from the first to the second line, we have used the fact 
that dR/dT k, — Stt close to the surface but at a finite distance below it, for a 
generic function f{R) = R — p.'^/R + eR?. To see this, one can solve (jH) and obtain 
R = — 47rr± (3^'' + IBtt^T^)^/^. Choosing the positive sign in order to have a positive 
cosmological constant in vacuum, one has p? — Rq/^/S ~ 10^'*'^. Then, even very 
close to the surface, one has \T\ ^ p^ and R w SttT. 
Integrating the Euler equation p' = —A'{p^ 

5/2 



hence 



Therefore, 




p)/2 just below the surface one gets 

5/2 



-3/2 



-3/2 



(19) 



(20) 



(21) 



V'^out - 2mtot V'^out 

To calculate the ratio given by pT|) . let us first consider the general case 
f{R) — R — p*/ R + eR^. We stress again that one generically expects the presence of 
the term ei?^, because there is no first principle from which to derive the functional 
form of f{R), and one has to think of it as the series expansion of an unknown f{R) 
coming from a consistent high-energy theory of gravity. As can easily be seen from 
([H]), the quadratic term does not influence the vacuum value Rg of the curvature 
scalar, which acts as the effective cosmological constant. Basically for this reason, 
the quadratic term is essentially unconstrained by cosmological data and solar system 
tests only allow weak constraints to be placed on it [17]. Taking now e ^ 0.1 in 
our units (a value several orders of magnitude smaller than the upper limit coming 
from solar system tests [17]), just below the surface we have dF/dR « 2e « 0.2 and 
Ci « — 2/Jo ~ —3/2 (because F' near to the surface). From ([^T|l . one then obtains 

-3/2 

2 V — _ ^ (22) 



^smg 



y/rout - 2mtot V^out - r 
In the case of a neutron star as modelled with a polytropic EOS (k 
and mtot ~ 1), one therefore has 



4, rout 



acR 



0.02(rout - r)-i/2 



10 



(23) 
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and the ratio |asing/aGRl is large at distances below the surface at which the fluid 



approximation is certainly valid. For instance, lasing/^GRl 20 for rout — r ^ 10 ^ 
1.5 mm, losing/aGRl ~ 600 for Tout - r 10"^ 1.5 /im, |asi„g/aGR| ~ 2 x 10^ for 



-^sing/ 

1" o 1 n I T 7- T ro MiH H/^T"mc^~nn/^j:i I/, . , 

3r| - 600 for rout - r ^ 10~^ 
rout — r ~ lO"""^^ ~ 1.5 nm. Note also that the ratio |asing/aGR| scales proportionally 
with the value of e, which we have taken, as already mentioned, to be several orders 
of magnitude smaller than the upper limits coming from solar system tests (17j . 

Let us now consider instead the case f{R) = R — n'^/R, as used in [16]. First, we 
need to evaluate a. Noting that ci — 2/Fo — —3/2 as r — > rout, because F' = at 
the surface [13], we have ci ~ 1 just below the surface. In order to see what happens 
instead at a finite distance below the surface, note first that dF/dR « —2Rq/{3R^). 
As already mentioned, solving ([6]) and imposing that the cosmological constant 
A = i?o/4 in vacuum is positive, one gets R = -inT + + len^T^y/^ « Snp 
for p 2> Rq- We can then write dF/dR « -~2/[3{8tt)^]RI/p^ for p > Rq, and therefore 
F' = {dF/dR) [dR/dT) {dT/dp) p' - p' {Rq/pY/Rq. Because Rq - 10-'*^ it is clear 
rF' <C F even at finite distances below the surface. From it then follows that 
ci w — 2/Fo = —3/2 also at finite distances below the surface. 

Let us now evaluate (|2T]) for p ^ Rq: using ci w —3/2 and dF/dR w 
-2/[3(87r)3]i?2/p3^ j|. becomes 

37r /5 



^smg 



acR 



25(87r)3 V2 



R^K'^m^^r out {rout -2mtot ) (^out -r) ^ . (24) 



It is therefore clear that tidal forces become increasingly more important, even in the 
particular case f{R) = R — p'^/R, for spheres with larger radius. As such, even for 
this particular form of f{R), the lengthscale on which the tidal forces diverge is much 
larger than the lengthscale on which the fluid approximation is valid, if one considers 
sufficiently diffuse systems: some examples are worked out in the Appendix. 

In conclusion, we have shown that the fluid approximation is still valid on the 
scale at which the tidal forces diverge just below the surface of a polytropic sphere 
in the case of the generic functions f{R) likely to arise in practice in a cosmological 
scenario. Even in the special case considered by Kainulainen et al. |16j , this continues 
to hold for configurations which are sufficiently diffuse. 



3. Physical and mathematical nature of the problem 

3.1. Differential structure and cumulativity 

It is clear from the above that the nature of the problem discussed here does not 
lie in the fluid approximation or in the specifics of the approach followed in |14| . 
but is related to intrinsic characteristics of Palatini f{R) gravity. These concern the 
differential structure of the action ([T]) and the resulting field equations. 

We recall that the Lagrangian of the action ([T]) is an algebraic function of 
R = g^^^Rf^i, and that R^^ is constructed from the independent connection T^^^. 
In more detail, 

R'^.ax = -dxT'^ua + d,T\^ + F'^^^F",, - F'^^^F",, , (25) 
and contracting the first and the third index one gets (see [H H] for further details) 

Rfiu ~ R^fj,\v — ^aT^i/ ~ ^i^r^A + r^Ar^iy ~ r'^'o-i^r^A- (26) 

It follows from this that f{R) has no a priori dependence on derivatives of the metric. 
Also, R depends only linearly on the first derivatives of the connection i.e., at least 
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in the case where / is hnear in R (which leads to GR), there are no dTdV terms 
(indices suppressed) as there would usually be in a field theory! One might expect 
that allowing / to be non-linear in R would introduce dVdV terms and solve this last 
problem, but we will see shortly that this is not the case. Note also that, since the 
metric has no a priori relation with the connection, one is dealing with a field theory 
with two independent fields, and so one cannot argue that having no quadratic terms 
in the connection is expected because the connection already includes derivatives of 
the metric. 

This lack of dynamics in the action is also mirrored in the field equations ^ and 
([3]). Variation with respect to the metric leads to ([2]), which includes no derivatives 
of the metric. As already mentioned, contraction of ([2]) gives ([6]), which algebraically 
relates R and T for a given f{R). Variation with respect to the connection leads to Q, 
after some integration by parts to "free" the connection. For a linear function f{R), 
this equation is just the definition of the Levi-Civita connection. When f{R) is non- 
linear, instead, ^ seems to include second derivatives of the connection. However, 
this is misleading because R can be completely eliminated in favour of T by using 
([6]) and, therefore, ([3]) can be trivially solved to give the connection as a function of 
the metric and the matter fields. As already mentioned, following these steps one 
can completely eliminate the connection in favour of the metric and the matter fields, 
and turn ^ and ^ into the single- field representation ([7]). This representation of 
the theory is more convenient and more familiar for discussing the dynamics. It also 
highlights once more that the metric fully describes the geometry, which is indeed 
pseudo-Riemannian, and that the independent connection is just an auxiliary field 

It is also interesting to note that one could introduce an auxiliary scalar (j) = F 
and re-write ([7]) as 



- ~ (^(V^0)(V.0) - ig^.(V0)2^ , (27) 

while, setting V(0) = i?0 — /, ® can be re-written as 

2V{(j)) ~ (j)V'{(j)) ^8nT. (28) 

Expressions (P7)) and ([^5]) are the field equations of a Brans-Dicke theory with Brans- 
Dicke parameter luq = —3/2, i.e. a theory described by the action 

^ ^ d^x^ (cf>R + T^a^^a'^^ - V{c^)] + Sm{9^.,^), (29) 



IGtt J ^ ^ y-r 20 

(see also [3 [26l [27] for more details about the equivalence of Palatini f{R) gravity 
and (jjQ = —3/2 Brans-Dicke theory). 

Returning to ([7|) , we note that this is a second order partial differential equation 
in the metric, just as in the case of GR, but that the left hand side includes up to 
second derivatives of F and consequently of T [i^ = FiR) and R — R{T)]. Usually, 
the matter action includes derivatives of the matter fields ijj (if the equation of motion 
of the matter fields is to be of second order, the matter action has to be quadratic in 
the first derivatives of the matter fields). Therefore, generically one has T = T{dtp, 
implying that (O includes up to third derivatives of the matter fields! 

In GR and in most of the proposed alternatives to it, the field equations include 
only first derivatives of the matter fields. The higher differential order in the metric 
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with respect to the differential order in the matter fields guarantees that the metric 
comes as an integral over the matter fields. Therefore, any discontinuities in the matter 
are "smoothed out" and are not inherited by the geometry (cumulativity of gravity) . 
We recall that in general the metric is not allowed to become a delta function or a step 
function (although the latter is allowed if no Dirac deltas are produced in the field 
equations, i.e. if the metric, in spite of being discontinuous, is a solution of the field 
equations in the sense of distributions: see for instance [28], section 3.7). However, 
this is clearly not true in Palatini f{R) gravity or in = — 3/2 Brans-Dicke theory, 
since the differential order of the field equations in the matter fields is actually higher 
than in the metric, implying that the latter is not necessarily an integral over the 
matter fields, but can be algebraically related to the matter fields and even to their 
derivatives. Because of this, a discontinuity in the matter fields or in their derivatives 
can lead to unacceptable singularities. A similar behaviour has been demonstrated 
in the post-Newtonian limit of the theory, where the post-Newtonian metric becomes 
algebraically dependent on the matter density fT3]. 

This unusual differential structure of Palatini f{R) gravity is at the root of the 
surface singularities discovered in [M] . The polytropic description of matter was used 
in [H] only because this made it possible to find analytic solutions and demonstrate 
the problem without resorting to numerical techniques. In fact, a more detailed 
description of the matter would make the problem even more acute. To see this, 
note that in the case of a perfect fluid one has T^^ = T^v{p,p), i.e. the stress-energy 
tensor does not include any derivatives, unlike the case of a microscopic description 
of matter. The fluid approximation actually "smoothes out" the matter distribution 
with respect to the microscopic description. This also explains why no singularities 
appear for 1 < F < 3/2: these values of F give a smooth passage from the interior to 
the exterior. In conclusion: abandoning the fluid approximation would just increase 
even further the differential order of the field equations in the matter fields and make 
it easier for singularities to appear. 

As a further confirmation that the introduction of microphysics cannot solve 
the problems caused by the algebraic dependence of i? on T [see ([5])] or, in the 
equivalent action ([55]), by the algebraic dependence of on T [see ([^ ]. let us note 
that this feature of Palatini /(i?) gravity introduces corrections to the standard model 
of particle physics already at the meV energy scale (see [M] and [29]). Both the 
calculation of [M] and that of [29] are performed in the Einstein frame. Although 
the use of the Einstein frame has been criticized [30j|l . this frame is equivalent to 
the Jordan frame and both are perfectly suitable for performing calculations |26j . 
However, one should remember that particles in vacuum follow geodesies of the Jordan 
frame metric, so this is the metric which becomes approximately Minkowski in the 
laboratory reference frame. This makes the Jordan frame calculation simpler and 
more transparent than the one in the Einstein frame. For this reason, and in order to 
highlight once again the problems caused by the algebraic dependence of R on T, we 
briefly redo the calculation of [26] and [29] in the Jordan frame. Let us first consider 
the equivalent action ([29]) and take the matter to be represented by a scalar field H 

X If the independent connection is allowed to enter the matter action, the results of I26| and I29| 
will of course cease to hold, as pointed out in [301 . In this case, also the surface singularities that we 
found in 1141 may disappear (see also section [3.21 . However, such a theory would be a generalization 
of Palatini f{R) gravity [see the action |(TJ], known in the literature as metric affine f{R) gravity [6]. 
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{e.g., the Higgs boson), the Lagrangian of which reads 
1 

2hV -^--"^ 



^[g^^d,Hd,H-'-^H'] (30) 



(we recall that we are using units in which G = c = 1). The vacuum of the action 
(EH) with ^ and f{R) =R- fj^/R [which implies ¥{((>) = 2^2(0 _ i)i/2] ^an easily 
be found to be 7f 0, (?i = 4/3 [the solution of ([28]) with T ^ 0] and 

g,.d.^d.^ ^-{l- - dt^ 

+ dr^/ (^l - _ -,_ r^dn^ « rj^^dx^'dx'' (31) 

(which is indistinguishable from the Minkowski metric for the purposes of a particle 
physics experiment because /i^ ^ A and r « ^Eanh)- One can then expand the action 
to second order around this vacuum (as usual the first order action is identically zero 
because the field equations are satisfied to zeroth order). However, it is easy to show 
that perturbing (|28p one gets Scf) ^ T/^"^ ^ m^SH^ / {h^ fi^) at energies lower than 
the Higgs mass (mn ~ 100 — 1000 GeV): replacing this expression in the action (P^)) 
perturbed to second order one immediately gets that the effective Lagrangian for the 
Higgs scalar is 



m^SH^ m\{d6Hf 



(32) 



at energies k ^ mn- At an energy fc = 10 eV (corresponding to a lengthscale 
L = h/k = 2x 10-4 m), using the fact that //^ _ A - 1/{H^^)'^ (where H^^ = 4000 
Mpc is the Hubble radius) and 6H mn (because k <C mn) and remembering 
that we are using units in which G = c = 1, it is easy to check that the first 
correction is of the order m'^SH'^ / {^^ hP) ^ {HQ^/\n)'^{mii/Mp)'^ 3> 1, where 
Ah = h/mn ^ 2 x 10^^^ — 2 x 10~^^ m is the Compton length of the Higgs and 
Mp ^ ^1/2 ^ (^cVG)i/2 ^ 1.2 X 10^^ GeV is the Planck massH Similarly, the second 
correction is of the order m^{dSH)^ / (p^h^) - (H^^ /Xn)^{mn/Mp)^{HQ^ /L)^ > 1. 

Note that replacing Scj) ~ m^5H^ / {'h? ji^) in ([29| gives also that the coupling of 
matter to gravity is described by the interaction Lagrangian 

m^SHW^ dHg\ mlSH\ ^ /^-i^^n 



1 



L 



(33) 



It is therefore clear that also the coupling to gravity becomes non-perturbative at 
microscopic scales. This is, once again, a consequence of the algebraic dependence of 
(j) on T, encoded in ((28|) , and this is in agreement with the singularities that we discuss 
in this paper. 

§ Equivalently, one can write the first correction as a self-interaction term m'^SH'^ / (p?h^): restoring 
the dependence on G this term becomes Gm'^&H'^ /(p?ti^). In "particle physics units" h = c= 1, the 
coupling constant is dimensionless and is given by Gm^/ ~ (^^q^"'^/Ah)^('tih/A^p)^ 2^ 1- 
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3.2. Overcoming the problem 

In section lXTj we have traced the root of the problem: it hes in the awkward differential 
structure of the field equations, in which the matter field derivatives are of higher 
order than the metric derivatives. This introduces non-cumulative effects and makes 
the metric extremely sensitive to the local characteristics of the matter. With this 
in mind, it is not difficult to propose a possible way out. Clearly, one would like 
to restore the cumulative nature of gravity. This requires the introduction of more 
dynamics into the gravitational sector of the theory. As an example of how to introduce 
more dynamics, let us consider a theory described by the action: 



where a should be chosen so as to have the correct dimensions. The term R^^^ R^i, is 
quadratic in the derivatives of the connection [see ([26]) ]. This implies that the action 
(|34p. even though linear in both R and R'^'^Rf^^, is quadratic in the derivatives of the 
connection and will not lead to an algebraic equation for the connection, unlike the 
action ([T]). Indeed, the field equations that one derives by varying the action ([51)) with 
respect to the metric and the connection are, respectively, 



and (|36p cannot be algebraically solved for the connection. Also, this theory cannot 
be re- written as an ojq = — 3/2 Brans-Dicke theory. In summary, a theory described 
by action does not seem to be sharing the unwanted characteristic of Palatini 
f{R) gravity: that after eliminating the connection, one ends up with the matter field 
derivatives being of higher order than those of the metric. In particular, R cannot 
be expressed as an algebraic function of T though the trace of eq. ([35)l . as in the 
case of eq. (I6|), nor can the independent connection be algebraically expressed simply 
in terms of the metric and derivatives of the matter fields (therefore introducing the 
higher differential order in the matter fields when it is replaced in the field equations). 

Theories with higher order invariants in the action, such as R^^R^'^ , have recently 
been considered in the Palatini formalism in [31]. Clearly, a more detailed analysis 
of the dynamics of such theories is needed in order to show in a clear way whether 
they exhibit the problem discussed here or other viability issues. Here we have used 
them solely to demonstrate that it might be possible to overcome the issues discussed 
here by generalizing the action. This clarifies the following point: such shortcomings 
are not generic to Palatini variation, but seem to be a specific problem of Palatini 
variation when used with the specific choice of /(i?) actions. 

4. Conclusions 

In this paper, we have discussed in detail the issue, raised in [14] . of surface curvature 
singularities appearing for polytropic spheres in Palatini /(i?) gravity. Simple 
gedanken experiments lead us to conclude that the presence of these singularities 
casts serious doubts on the viability of the gravity theory. Concerning the objection, 
raised in [16], that polytropic EOS's are too idealized to allow one to rule out Palatini 
f{R) gravity, we stress that among the EOS's not giving a regular static spherically 




(34) 




(35) 
(36) 
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symmetric solution [14] there are perfectly physical cases such as a degenerate non- 
relativistic electron gas or an isentropic monatomic gas. Regular solutions for these 
configurations exist even in Newtonian mechanics, and we have argued that a theory 
not providing such solutions should be considered, at best, as being incomplete and 
as being disfavoured for giving viable alternatives to GR. We have also presented 
quantitative results for the magnitude of the tidal forces exerted just below the surface 
of polytropic spheres, showing that, for generic forms of f{R), the lengthscale on 
which the tidal forces diverge due to the curvature singularities is much larger than 
the lengthscale at which the fluid approximation breaks down. This generalizes the 
calculation of Kainulaincn et al. [16]: while their result (that the tidal forces only 
diverge at lengthscales on which the fluid approximation is not valid) is correct for 
their particular choice of f{R) and for neutron stars, we find that it does not hold, 
even with their choice of f{R), for more diffuse configurations and does not hold, 
even in the case of neutron stars, for more plausible choices of f{R)- Finally, an 
analysis of the differential structure of the field equations for the theory has been 
presented, which sheds light on the origin of the problem, showing that the appearance 
of singularities is not related to the fluid approximation. On the contrary, abandoning 
the fluid approximation would make the problem even more acute. The addition of 
more dynamics to the theory seems to be a potential way out of this difficulty. 

Appendix 

In this appendix we evaluate - which gives the ratio Osing/flGR in the special case 
of f{R) = R — iJ-'^/R - in several contexts. Our calculations will show that even with 
this special choice of f{R), the fluid approximation is still valid on the scale at which 
the tidal forces diverge if the configuration under consideration is sufficiently diffuse. 
As in section [2^ we use units in which G = c = Mq = 1. 

Let us first consider a solar type star with mass rristar ~ 1 surrounded by 
a gas cloud with mass rricioud ~ 0.01 and radius rout ~ 10^** km, composed of 
monatomic isentropic gas. This is a perfectly physical configuration, although possibly 
not an astrophysically fully realistic one (note however that rout ~ 10^^ km is 
approximately the outer radius of the Oort cloud). The total mass of such a system 
is mtot ~ 'Tigtar + ™cioud ~ Ij ^^d the polytropic constant of the cloud is k w 9 x 10^^. 
From (pi)) . one gets 

«6x 10™(rout-r)-^ (37) 

and the tidal force becomes comparable to that of GR at a distance below the surface 
comparable to rout! Taking, for instance, a value of rout — r « rout/10 ~ 10^'^, the 
tidal forces in Palatini f{R) gravity would be 6 orders of magnitude larger than in 
GR. At this distance from the surface, the mean distance between the particles of the 
fluid is i 1/n^/^ « {mp/p)^/^ « 10~^ (where n — p/mp is the number density and 
nip ^ 10^^^ is the mass of the proton). An upper limit for the MFP at this distance 
from the surface can be calculated assuming a cross section a ^ (1 A)^ « 5 x 10~^^, 
giving ^MFP 1/('T-^) ™p/(po') ~ 10^"^ ^ ''out — ''[jjl- Also, note that the average 

II This is an upper limit because it assumes a "geometrical" cross section for encounters between the 
atoms (1 A is approximately the size of a hydrogen atom). However, for hydrogen- hydrogen collisions 
in the lab, cr(HH) ~ 20 A'^ |20] while, for instance, hydrogen-lithium collisions have cross sections 
which are about 1200 |21| . More importantly, if the fluid is (even partly) ionized, the cross section 



Rising 

acR 
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velocity of the particles in the cloud can be evaluated from p — nk^T (k^ and T 
being Boltzmann's constant and the temperature) using and (PO)) . and is Uav ^ 
{mtot{rout-r)/[rout{rout-'^intot)]V^'^- For rout-?' ~ ?'out/10, one has Uav 4x 10"^, 
which is comparable with the virial velocity Wviriai ^ (?7itot/r)^^^ ~ 10~^, and so the 
polytropic coefRcient k needed to support the cloud is plausible. In conclusion: for 
this configuration, the lengthscale on which the tidal forces in Palatini f{R) gravity 
are larger than in GR is certainly larger than the lengthscale on which the fluid 
approximation is valid, whether this scale is taken to be the mean interparticle distance 
or the MFP. 

We would argue, however, that the relevant scale is actually the interparticle 
distance, because we are considering here static equilibrium configurations. A way to 
understand this point is to consider how one derives the hydrodynamic equations from 
the Vlasov equation (i.e., the conservation equation for the phase-space distribution 
f{x,v) in the case of a collisionless fluid, that is, one with infinite MFP; see for 
instance P3], paragraph 27): 

^ + Ky^ + ^E1=Q (38) 

where v'^ = dx^ /dt is the velocity, m is the mass of the particles and is the force 
(thought of as dependent only on position and not on velocity). By integrating over 
all velocities, one easily obtains the mass conservation equation dtp + V • [pv) — 0, 
where p — m J fd^v is the density and = m J fv'^d^v/ p is the average {i.e., 
macroscopic) velocity. Similarly, one can multiply pSJl by and integrate over 
all velocities. If the velocity distribution is isotropic 131 one then gets the Euler 
equation {dt+v ■'V)v = — (Vp)/p + F/m, where one uses the isotropy of the velocity 
distribution to define the pressure as pSij = rn J /(u* — v^){v^ — v^)d^v. If the phase- 
space distribution / is specified, the mass conservation and the Euler equations are 
clearly a closed system of equations. Therefore, one does not need to consider higher 
order moments of the Vlasov equation, and the system under consideration is a fluid 
in spite of the MFP being inflnite. A typical example of this situation is, for instance, 
that of Dark Matter in a Friedmann-Robertson- Walker universe. Since the velocity 
distribution is isotropic because of the cosmological principle. Dark Matter can be 
treated, at the background level, as being a fluid (c/. for instance (3.10)-(3.12) of 
|24| ) . Similarly, for a gas trapped in a box, the fluid approximation is valid on scales 
larger than the interparticle distance, whereas the MFP can be infinite if the fluid is 
non-coUisional. This is indeed the case considered in the textbook derivation of the 
perfect gas law (see, for instance. Landau, Lifshitz and Pitaevskii [25], chapter 4), 
where the only necessary hypothesis is the isotropy of the velocity distribution. The 
role of the box is, in our case, played by the gravitational potential well. 

can be much larger, because Coulomb forces are long range (in strongly coupled plasmas it is actually 
common to have a MFP shorter than the interparticle distance |22|). 

^ Of course, one may object that if the fluid is collisionless there is no interaction which can make 
the velocity distribution isotropic. However, we are interested here in showing that for equilibrium 
configurations the MFP has nothing to do, from the conceptual point of view, with the lengthscale 
at which the fluid approximation breaks down. Moreover, one can always think of a tiny interaction 
between the particles (resulting in a huge MFP) which can make the velocity distribution isotropic 
in a sufficiently long time (comparable with the mean free time). From the conceptual point of view, 
one can also think of shooting a beam of collisionless particles into a box (or a potential well): the 
initially focused velocity distribution will become isotropic due to the small irregularities in the walls 
of the box or in the gravitational fleld. 
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Let us now consider a polytropic sphere with rritot ~ 0.1 and rout ' 
(J?0 being the radius of the Sun). The polytropic constant is then k 
becomes 



lOHrout-r)- 



2OOi?0 « 10^ 
i 2 X 10"^, and 



(39) 



At a distance of ~ 1.5 km below the surface, therefore, tidal forces are ~ 10^ times 
stronger in Palatini gravity than in GR, while the forces in the two cases become 
comparable at a distance > 15 km. Now, from (PO)) . at a distance rout — T ~ 1.5 km « 1 
we have p w lO^'^^ ~ 10®i?o- Although Tout — r is certainly smaller than the upper 
MFP limit introduced above, at this density the mean distance between the particles 
of the fluid is £ ^ l/n^/^ ~ (nip/p)'^/^ « 10"'^ ~ 0.1 mm. 

The same considerations apply, although marginally, for a polytropic sphere with 
?Titot ~ 1 and Tout ~ ^0 ~ 5 X 10^. The polytropic constant is then k w 2 x 10^, and 
(l24l) becomes 



^smg 



I0-^%rout - r)' 



(40) 



from which it follows that the difference between the tidal forces becomes important 



for Tout — r < 7 X 10 



1 cm. At this distance below the surface the density 



is /9 ~ 5 X 10 ^'^ ^ 10^^i?o, and the mean distance between the fluid particles is 
e « l/ni/3 ~ (TOp/p)i/3 10-8 _ 0.01 mm. 
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